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Abstract 

In this work we consider a coupled system of Schwinger-Dyson equations for self-energy 
and vertex functions in QED^. Using the concept of a semi- amputated vertex function, we 
manage to decouple the vertex equation and transform it in the infrared into a non-linear 
differential equation of Emden-Fowler type. Its solution suggests the following picture: 
in the absence of infrared cut-offs there is only a trivial infrared fixed-point structure in 
the theory. However, the presence of masses, for either fermions or photons, changes the 
situation drastically, leading to a mass-dependent non-trivial infrared fixed point. In this 
picture a dynamical mass for the fermions is found to be generated consistently. The 
non-linearity of the equations gives rise to highly non-trivial constraints among the mass 
and effective ('running') gauge coupling, which impose lower and upper bounds on the 
latter for dynamical mass generation to occur. Possible implications of this to the theory 
of high-temperature superconductivity are briefly discussed. 
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1 Introduction 



Three-dimensional Quantum Electrodynamics {QED3), with an even number of 
fermion flavours, apart from serving as a toy model for studying chiral-symmetry 
breaking and confinement, also constitutes a physically interesting theory per se, 
in view of its possible applications in modelling novel (high-temperature) supercon- 
ductors 

Chiral symmetry breaking or, equivalently, dynamical mass generation for fermions, 
in even-fiavour QED3 has still many unresolved issues. One of those is the exis- 
tence of a (dimensionless) critical coupling, above which dynamical mass generation 
for the fermions occurs In the context of large- iV treatment, which at present 
constitutes the only well-studied approach, the role of the dimensionless coupling [] 
is played by the inverse of the fermion fiavour number N. The issue of the exis- 
tence of a critical coupling in QED^ is a delicate one 0. Many of the original 
approximations leading to its existence have been questioned, in particular the 
fact that wave-function renormalization effects have not been properly accounted 
for. Recently, however, the incorporation of such effects, still within the large -A^ 
context, appears to corroborate [0]-0 the qualitative picture advocated in In 



addition, the latter is also supported by lattice simulations [|T0 . 

Nonetheless, the situation is far from being conclusive. The fact that the critical 
N, below which dynamical mass generation occurs, is found to be of order 3 (in a 
four-component notation for the fermions) provides motivation for searches beyond 
the large- treatment. Moreover, up to now the gauge coupling in the infrared has 
been treated as an arbitrary parameter, whose size has not been restricted by an ad- 
ditional dynamical constraint. At present we are lacking a self-consistent treatment 
of the dynamical Schwinger-Dyson (SD) equations involving the vertex function on 
an equal footing with the self-energy and gap functions. In all the approaches so far, 
at least within the large A^ treatment that we are aware of, one invokes a specific 
Ansatz for the vertex, by the sole requirement of satisfying some truncated form 



of the Ward-Takahashi identity stemming from gauge invariance 0-111]. The lack 
of dynamical information for the coupling poses problems; for instance, its size in 
the infrared is treated as an arbitrary parameter, being assumed to merely exceed 
a critical value, if one wishes to trigger chiral symmetry breaking. 

Another point, related to the above, which is already familiar from studies in 
the case of four-dimensional non-Abelian gauge theories, is whether chiral-symmetry 
breaking is associated with confinement of charges |T^. This issue acquires physical 
importance in view of the condensed-matter applications. In particular, it may shed 
more light in the dynamics of spin-charge separation, by analogy with the physics 
of strong interactions [|13[ . 



In this work we shall not deal with issues of confinement, which exists in QED^ 



^In three dimensions the couphng has dimensions of mass. One can still define, however, 
dimensionless couplings by dividing with a dynamically generated scale, which in the large iV- 
treatment arises by demanding that ||] e^N = 8a, with the scale a kept fixed as — ^ 00. The 
dimensionless coupling is then defined as the ratio e'^/8a = 1/N. 



despite its abelian nature. Instead, we shall attempt a novel approach to chiral 
symmetry breaking, independently of a large N treatment, by studying the coupled 
fermion and photon self-energies and vertex SD equations in the context of a method 
first introduced for the case of four- dimensional QCD The novel ingredient is 
that we concentrate on the semi-amputated vertex, defined in section 2, which is 
the correct gauge-invariant quantity to determine a physically meaningful "running" 
coupling ('effective charge'). The fact that QED^ is superrenormalizable in the 
ultraviolet does not preclude the possibility of defining such a quantity, having non- 
trivial structure in the infrared. As we show in sections 2 and 3, the existence of a 
non-trivial infrared fixed point is a property only of the infrared-regularized theory. 



as conjectured in [0. In other words, in the absence of fermion (or photon) masses, 
the infrared singularities of the vertex equation will force the effective charge to 
vanish at zero momentum transfer, thereby excluding the possibility of an interesting 
infrared behaviour. From a condensed- matter application point of view this would 
correspond to what is usually called the Landau- Fermi liquid theory (trivial infrared 



fixed point) |TT1. In the presence of masses, and in particular fermion masses, we 



show in section 3 that there is a non-trivial infrared fixed point structure, stemming 
from the fact that the effective charge obtained as a self-consistent solution of the 
n on- linear vertex SD equation, is driven to a finite positive value, which can be large 
enough to trigger dynamical generation of a fermion mass. This implies that the 
phenomenon of chiral symmetry breaking is intimately associated with deviations 
from the trivial infrared fixed-point structure. 

We should stress that, as a result of the non-linearity of the vertex equation, there 
are delicate constraints between the fermion mass and the effective charge, which are 
responsible for the appearance of regions of the latter for which dynamical genera- 
tions occurs. At present, these restrictions appear as a consequence of mathematical 
self-consistency of the truncated equations. It is not clear to us, whether the up- 
per bounds on the effective charge, imposed by the present cubic approximations 
for the vertex corrections, will survive the inclusion of higher orders. In contrast, 
we believe that the lower bounds will survive such a treatment, thereby indicating 
the existence of a critical coupling above which dynamical mass generation will oc- 
cur. This is physically appealing, given that one would not expect a weakly-coupled 
theory to be capable of breaking dynamically chiral symmetry. 

The layout of this article is as follows: in section 2 we set up the SD equa- 
tions that we wish to study, and discuss in detail the approximations employed. 
We demonstrate that, under certain assumptions to be justified retrospectively in 
section 4, the equation for the semi-amputated vertex decouples from the rest, and 
hence can be solved separately. Moreover, we establish the absence of a non-trivial 
infrared fixed point rigorously (within the cubic approximation for the vertex), by 
casting the SD equation for the effective charge in a form of a non-linear differ- 
ential equation, known as Emden- Fowler equation |TB, |r7|. In section 3 we study 



the equation for the vertex in the presence of a fermion mass, acting as an infrared 
regulator. We derive the appropriate non-linear differential equation describing the 
infrared behaviour of the running coupling, and solve it to demonstrate the ex- 



istence of a non-trivial infrared fixed point. The so-derived running charge is a 
monotonically decreasing function of the momentum, tending asymptotically to a 
constant positive value in the ultraviolet. In section 4, we solve the equations for 
the photon and fermion self-energies, which in our approximations decouple from 
each other and depend only on the semi-amputated vertex function. We then verify 
the self-consistency of the approach. In section 5 we examine the self-consistency of 
the dynamical generation of the fermion mass, by solving the appropriate SD equa- 
tion upon substituting the solution for the semi-amputated vertex found in previous 
sections. The self-consistency of the approach restricts the allowed regions of the 
effective charge, implying the existence of a lower bound (critical coupling) but also 
of an upper one. In section 6, we examine an alternative type of infrared cut-off, 
namely that of a (bare) covariant photon mass term. This case also exhibits a non- 
trivial infrared fixed-point structure but, in contrast to the monotonic decrease of 
the effective charge in the case of fermion masses, here the coupling initially increases 
in the infrared, then displays a local maximum, and eventually decreases, tending 
asymptotically to a constant value in the ultraviolet. Some possible applications of 
this behaviour, inspired by condensed-matter physics, are briefly discussed. Finally, 
in section 7 we present our conclusions and outlook. 



2 The SD equation for the semi-amputated ver- 
tex 

In this section we will first set up the SD equations for the photon and electron self- 
energies, and the photon-electron vertex; then we will define the semi-amputated 
vertex and derive its corresponding SD equation. As we will explain, the latter 
governs the behaviour of the effective coupling in the infra-red. The derivation of 
the SD equations for the photon propagator A^j,, the electron propagator Sp, and 



the photon-electron vertex proceeds following standard methods |T^, |1^ (see 
figures 

The full photon propagator A^j,, its inverse A~J, and the full vacuum polarisation 
Uf^i, in Euclidean space are related by 

= Aoi(?) + n^.(g) 

^oi(?) = Q^^f^u - (1 - ■^)qt^qu 
^tiu{q) = {q^ V - q^,qy)Ii{q) 

= is^u-q^qu/qW + q'^iq)]~' + Cq,qu/q' (2.1) 

where ^ is the gauge fixing parameter (in covariant gauges). The corresponding SD 
equation reads 

A-J{q) = Ao-i(g) + e'l ^^Tr[T,SpTA] + ■ • • (2.2) 



The SD equation for the electron propagator Sp is given by 



i2ny 

Finally, the SD equation for the photon-electron vertex has the form 



(2.3) 



7m -e 
+ . . . 



r.^F(fc-pi)F^^p(fc)F^A°^(fc+P2) 



(2.4) 



with pi + P2 + P3 = 0. The ellipses on the right-hand side of ( p.2|) -( p^ , denote 
the infinite set of terms containing the two-particle irreducible four-point function 



Ig, . Although we are not working in the context of a large analysis, we note 



that the above truncation is compatible with working to leading order in resummed 
expansion. 



+ 



Figure 1: Schematic form of the SD equation for the gauge field propagator in 
resummed perturbation theory. The blobs on the right-hand-side indicate the full 
(non-perturbative) vertex corrections. 




Figure 2: Schematic form of the SD equation for the full fermion propagator. Blobs 
indicate full non-perturbative quantities. Notice the full vertex appears on both 
ends of the internal photon line. 



We next define the scalar quantities A, B and Q as follows: 

1 



A{k) 



B{k)k'^ 



(2.5) 
(2.6) 



Figure 3: The SD equation for the vertex 



and 

r^(Pl,P2,P3) = Q{PuP2,P3hl^ ■ (2-7) 
The quantity B is related to 11 defined in (|2.1|) by B{q) = 1 + n(g). The definition 
in ( |2.(j| ) imphes that the longitudinal pieces of the photon propagator will be dis- 
carded in what follows. Of course, there are no rigorous field-theoretic arguments 
justifying their omission or inclusion. The correct treatment of such terms neces- 
sitates a formalism which would allow for the self-consistent truncation of the SD 
series in a manifestly gauge-invariant way; unfortunately, no such formalism exists 
to date. The standard lore when writing down SD equations is to use in ( |2.3| ) the 
form of A^^ given in (|2.1|) , setting ^ = (Landau gauge). While in QED^ this 
choice renders the vertex corrections unimportant in the ultraviolet, it appears to 
be less compelling in the context of the superrenormalizable QED^. In addition, 
it is known that, while the conventionally defined fermion self-energy and photon- 
fermion vertex depend explicitly on the gauge-fixing parameter ^, it is possible to 
construct -at least at one-loop - a ^-independent fermion self-energy and vertex, by 
resorting to the diagrammatic rearrangement of the S'-matrix known as the pinch 
technique It turns out that the ^-independent fermion self-energy and vertex 

so constructed coincide with their conventional counterparts, if we choose for the 
latter the special value ^ = 1 (Feynmann-t' Hooft gauge) [^. Furthermore, as has 



been formally shown in ||2^, all longitudinal pieces appearing in ( |2.1| ) vanish from 
physical observables, such as S'-matrix elements, to all orders in perturbation theory. 
Thus, one is led to a generalized form of the Feynman-t' Hooft gauge, known as the 
"stagnant gauge", where only the 6^i, part of the vacuum polarization contributes, 
to all orders in perturbation theory. This gauge will be adopted throughout the 
present article. 

Following [|l^ and [0 we define the semi-amputated vertex G as 

G{pi,p2,P3) = Z{pi,p2,P3)GiPl,P2,P3) (2.8) 

with 

Zip,,p2,P3) = B-'/\p,)A~'/\p2)A-''\p^) (2.9) 

This definition proves very useful in reducing the complexity of the set of equa- 
tions ( |2.3| )-(p^, under certain approximations to be discussed in detail in the next 
sections. In addition, the quantity 

9RiPl,P2,P3) = eG{pi,p2,P3) (2.10) 



provides a natural generalization of the concept of the running or effective charge in 
the context of superrenormalizable gauge theories [^], such as QED^. This running 
of the coupling should be understood as a Wilsonian rather than Gell-Mann-Low 
type, in the sense that it is not associated with ultraviolet infinities; instead, it 
expresses a non-trivial infrared structure of the theory [|l^ . 



Note that in four-dimensional QED {QED4) the effective charge e^g is defined 
in terms of the photon vacuum polarisation as: 

els{<l')=e'[l + Il{q')]-\ (2.11) 



and is a gauge-, scale-, and scheme-independent quantity |2^ to all orders in pertur- 
bation theory, e^g depends explicitly on and the masses of the fermions inside the 
vacuum polarisation loop. In the limit where the fermion masses can be neglected, 
egg(g^) coincides with the running coupling obtained by the /? function of QED4, i.e. 
the solution of the usual renormalization-group differential equation. An advantage 
of the definition given in (|2.1CI| ) is that it captures the running coupling even in the 



case of scalar theories [19|, where, due to the absence of Ward-Takahashi identities, 
the role of the gauge boson self-energies is not as prominent as in gauge theories. 
As we shall see in section 5, the interpretation of qr defined in ( p.lO| ) as a running 
coupling is also justified by the form of the SD for the fermion mass gap. 

The equation for the semi-amputated vertex G may be obtained from (|2.4| ) by 
multiplying both sides by the factor Z{pi,p2,P3), i.e. 

r d?k 1 1 1 

G(pi,P„P3)7. = ^(Pi,P2,P3)7. -e'j (^^V^7.^7.^^^ (2.12) 



where 



= G(p3, k + p2,Pi - k)G{pi, -k, k - pi)G{p2, k, -k - P2) (2.13) 



In what follows we shall restrict ourselves to the case where the photon momentum 
is vanishingly small, and thus one is left with a single momentum scale p. One can 
then define a renormalization-group f3 function from this "running" coupling G{p) 
by setting 

P^p±G{p) (2.14) 
dp 

In order to further simplify the SD equation for G{p) we make the additional approx- 
imation that & = G^{k), i.e. a cubic power of a single G{k) depending only on the 
integration variable k. This approximation will be justified by the self-consistency 
of the solutions. 

Carrying out the gamma-matrix algebra using the formulae 7;_t7^ = —d, and 
Ifilplfi = {d — 2)7p valid for 4 x 4 gamma matrices in d{= 3)-dimensional Euclidean 
space, one obtains in a straightforward manner: 

1 r d^k s 1 1 



Several remarks are now in order. First, one observes that Z{p) 1 for p ^ oo, 
where perturbation theory is vahd. This is inferred from the fact that in such a case, 
as can be readily verified, the functions A{p),B{p) —^1 + 0{e^/p). In addition, in 
the ultraviolet region, p — oo, gauge invariance requires Q{j)) ~ Second, 
from (|2.15|) one observes that, if G stays positive, which is expected for any physical 
theory, then, as a result of the positivity of the integrand, G{p) > Z{p) for any p. 
Thus, one has the following basic properties of which stem directly from the 

integral equation ( |2.15|) : 



G{j>) > Z{p), for all p, ; G{p) ~ B'^^^p) ^1, p ^ oo (2.16) 

^ 1/2 

Notice that G{p) in the ultraviolet is thereby given by (1 + n(p)) as in QED^. 
This is the perturbative result, which, as we shall see later, is modified non-trivially 
in the infrared, in a way consistent with gauge invariance. Moreover, as we shall 
see later, the self-consistency of the approximations employed will require G(0) >> 



3/2. 

Our next assumption is that in the infrared regime, k/a « 1, which is of 
interest to us here, the effects of the inhomogeneous term Z{p) can be ignored. This 
assumption will be justified later on, when we consider a non-trivial self-consistency 
check of the solutions. We now remark that by ignoring the effects of Z{p) one can 
decouple the equation for the (gauge-invariant) amputated vertex from the equations 
for A{p),B{p). As we shall discuss in subsequent sections, these latter equations also 
decouple from each other, depending only on the vertex function G{p). 

Because of this, we commence our analysis from the SD equation for the vertex 
G(p), which we solve upon ignoring the effects of the inhomogeneous Z{p) term. 
Thus we arrive at the homogeneous equation 

This integral equation involves only one unknown function, namely G, which must be 
self-consistently determined. Note that this equation is invariant under the rescaling 
G G/e. This indicates a straightforward extension of the analysis to a large- 
treatment, given that N can be absorbed in a redefinition of e^. 

It is easy to see that, written in the form ( |2.17| ), the equation does not admit 



physically acceptable solutions, i.e. solutions with G > aiid finite^ Indeed, setting 
p = one obtains after the (trivial) angular integration 0: 



^Solutions that blow up in any point of the integration region are discarded. 

•^These arguments remain unaffected even in the presence of an inhomogeneous term Z{p), such 
that Z{0) is non-negative and finite. The non- negative nature of Z{p) ( p. 9 ) stems from that of 
A{p), which is guaranteed from general renormalization-group arguments 125 ] . 



Finiteness of G'(O) requires that the integrand of the right hand side of ( p.l8| ) con- 
verges at 1/ — i> and oo. The ultaviolet hmit does not present a problem, because 
the kernel vanishes like y"^, which is consistent with the superrenormalizability of 
the theory as well as the fact that the amputated vertex tends to 1. In the infrared 
limit y — > 0, however, the kernel blows up. For the integral to remain finite at that 
point, as required by the finiteness assumption for G(0), G^{y) must approach zero 
as y", a > 1/3, thereby implying that G{0) = 0. However for that to happen the 
integrand in ( p.l8|) must change sign, which would in turn imply that G{y) itself 
must change sign somewhere in y. According to our assumption above this is not a 
physically acceptable situation. 

To show rigorously that there are no physically acceptable solutions leading to a 
non-trivial infrared structure we next convert the integral equation into a non-linear 
differential equation of the type known in the mathematical literature as Emden- 
Fow/er equation ||T^, |T^. To this end, we perform the angular integration in ( ^^.17] ), 
to arrive at the equation: 

2 a dk ^^^ ^ 1^ + P| 



G{p) = —- G{ky\n\—^\ (2.19) 

37r"^ p Jo k k — p 

where we have set = 8a to make contact with the usual large- definition 0. For 
us, however, the number of fermion flavours is not assumed to be necessarily large. 
In fact, for brevity we set = 1 (in a four- component notation for the fermions) 
throughout this work. Next we introduce the dimensionless variables x = p/a 
and y = k/a. Since we are interested in the infrared behaviour of the model we 
consider the limit a; << 1, for which one obtains by expanding the logarithms in 
the integrand: 

A/ N 2 dy ^o, ,y + X, 

G{x) = / JLG\y)\n\'- 1 ~ 

irc^x Jo y y — X 

4 r , 4 f^dy^. 



Differentiating appropriately with respect to x, we arrive at the following differential 
equation for small x: 

x'^ + 3x'^ + ^^G%x) = x«l (2.21) 
dx"^ dx 37r"^ 

It is convenient to rescale G by setting: 



G.J^G (2.22) 



Then, eq. ( p.21| ) becomes: 



x^'L^ + ^:,^'^ + G\x) = ^, x«l (2.23) 
dx-^ dx 



Upon the change of variables ^ 



1 

2x^ 



G = 2'^/^ ^(0) the equation becomes of 



Emden- Fowler type |16, 17 



(2.24) 



As discussed in the mathematical literature ||T6|, the only non-trivial real solutions 
of ( p.24|) , as ^ — +00, are oscillatory about zero of simple sinusoidal form, which 
however oscillate infinitely rapidly as x — 0. The amplitude of the above solutions 
behaves like x"^/^, for x « 1. We interprete this behaviour as indicating an 
instability of the massless-fermion ground state. 

It is interesting to notice that the only power-law solution of ( p.23|) for x << 1, 
is purely imaginary, i.e. 



G{x) 



.V5 

^—x 



1/2 



X « 1 



(2.25) 



This solution would imply a 'trivial infrared fixed point structure' given that its 
associated (3 function vanishes at x = 0. However, the fact that ( p.25| ) is purely 
imaginary would again suggest instability. 

The above analysis constitutes a rigorous proof that, within the context of proper 
(i.e. finite and with finite-derivatives) solutions, and modulo the approximations 
discussed, no non-trivial infrared-fixed point is possible in QED^ in the absence of 
an infrared cut-off. This was conjectured in ref. [jl^], but here we have given an 
analytic proof. This motivates one to look for the existence of a possible non-trivial 
infrared fixed-point structure in the presence of fermion and/or photon masses. In 
the next section we shall discuss the case when the fermions develop a mass m. As 
we shall show, the existence of a non-trivial infrared fixed point is guaranteed due 
to the form of the resulting equations. 



3 Equation for the Vertex in the case of non-zero 
fermion mass 

As a first kind of infrared cut-off in the integral equation (|2.17|) we shall consider 
the case of a fermion mass gap m{p) = where S(j9) is the fermion self energy. 
In that case the fermion propagator Sp becomes: 

^-'*> = A(k)(lfl,„,(k)) 

For our purposes below we assume that m/(p) ~ m/(0) = 7^ 0. In that case the 
integral equation (|2.17|) becomes: 



Gip) = 77^ / 77;-^^G''ik)J—-—J—-—^+—^ / j^G^k) 



3 7 (27r)3 ^ \k^ + m}){k-py 3 7(27r)3 ^ \k^ + m})^ {k - p^ 

(3.2) 



Notice that the effects of the fermion mass are not given simply by just adding 
a mass squared term in the fermion denominators, but they result in additional 
structures in the integral equation. 

Performing the angular integrations one arrives at: 



dyf{y)ln 



y + X 



y-x 



G\y) 



where x = p/ a, m = vtijja are dimensionless, and 

y'^ + 3m^ 



!{y) = y 



(y^ + m 



2\2 



> 



Differentiation with respect to x yields: 



x—G(x) 
ax 



Stt'^x 



dyf{y) In 



y + X 



y 



X 



+ 



2xy 



x^ 



y' 



G\y) 



(3.3) 



(3.4) 



(3.5) 



One observes that formally as x — the right-hand-side vanishes, provided that 
G is finite. This indicates the existence of a fixed point. As we shall show below 
this is confirmed analytically by converting the integral equation into a non-linear 
differential equation. 

An additional feature which one would have hoped to study already at the level 
of the integral equation ( |3.5| ) is the monotonicity of G. Unfortunately the kernel 
in ( |3.5| ) is not manifestly positive to allow for such an analytic proof at the level 
of the integral equation for generic values of x, and one has to resort to numerical 
treatments, which fall beyond the scope of this article. However one can already 
infer from ( 3.5|) that, for high momenta x » 1, a monotonically decreasing G 
is consistent with the expectation that in this regime G is essentially given by its 
perturbative expression which asymptotes to 1. The analysis is omitted because it 
is straightforward. 

For low momenta, on the other hand, the behaviour of G{x) will also be shown 
to be monotonically decreasing, starting from a non trivial fixed point. This will be 
achieved by converting the integral equation into a differential one. Unfortunately, at 
present, we cannot analytically derive the monotonicity for intermediate momenta. 

To derive the differential equation from ( |3.3| ) we follow a similar analysis to the 
one leading to (^^ ). First, one expands the logarithms for small x << 1, thereby 
writing the equation as: 



G{x) 



5m 



37r2x^ Jo 
dyG{yf 



dyGivf 



yl _|_ jy^2 



dy 



+ 



{y^ + im?Y Stt^ Jx {y"^ + im?y 



y2 _|_ ^2 

dy A3 



G\y) + 



G\y) 



(3.6) 



Differentiating with respect to x one arrives at: 



x{x^ + m^f^G{x) + 3(x2 + m'f^Gix) + -^(x^ + 3m^)G\x) = (3.7) 
ax"' ax Stt^ 



The above equation can be solved numerically, to which we will come later on. 
However, in the infrared region x « m the equation accepts an analytic treatment, 
as we discuss below. In this region the equation ( P?7| ) is approximated by: 

x^G{x) + ^^G{x) + ^^G\x) = Q (3.8) 



It is immediate to see that a special power-law solution is given by (for positive 
G{x)) by: 

G{x) = rriTT^x-^/^ (3.9) 

Notice the infrared divergence of this type of solutions even in the presence of a 
(bare) fermion mass. The associated renormalization-group /3 function ( [2.14] ) for 
this case reads: 

f3{x) = --G ~ x"^/^ +00, as X ^ (3.10) 

indicating the absence of an infrared fixed point. The associated operator appears 
to be relevant (negative scaling dimension), which implies the possibility that the 
theory be driven to a non-trivial fixed point. 

However, in the infrared regime x << 1, one can find a different type of solu- 
tion [0: 

G = m7r^ — x^O (3.11) 
2v^l + c2x' ^ ^ 

where c is a constant of integration to be fixed by the boundary condition at x = 
implied by the integral equation, to be discussed later on. For physical solutions c 
is assumed positive. 

This type of solutions has a renormalization-group /3-function ( p.l4|) of the form: 



/? = -G{x) + -^^G\x) ^ - \ 0, x-0 (3.12) 

V37rmc (1 + c^x) 



from which we observe the existence of a non-trivial (non-perturbative) infrared 
fixed point at G* = ""^^^ > 0. Such a fixed point is the result of the dynamical 
generation of a parity- invariant, chiral-symmetry breaking fermion mass [^, indi- 
cating the connection of the phenomenon of chiral symmetry breaking in QED^ to 
a non-trivial infrared fixed point structure, in agreement with the expectations of 



ref. 11 



The non-trivial fixed-point solution ( p.ll| ) is not compatible with the integral 
equation ( p.2| ) for any value of the fermion mass m. Indeed one can derive a boundary 
condition for G{0) from (|3.2| ), which reads: 



^(0) = / dy^—^G\y) + / dy G\y) (3.13) 



In contrast to the massless case ( |2.18| ), G(0) is now a finite constant, Trmcy 3/8, 
as seen from and tliis allows for a compatibility of the solution (|3.11|) with 



T3D , provided that mG'(O) satisfies certain conditions to be specified below. 
To this end, we split the y-integration in ( |3.13|) into two intervals: (i) y G [0, m) 



where the form ( p . 1 1[ ) is valid to a good approximation, and (ii) y G [m, oo), where 
we approximate G{x) by its perturbative asymptotic value G ~ 1. By this latter 
approximation we overestimate the actual value of the integration, given that G is 
actually slightly smaller than unity for finite p, approaching it only asymptotically 
(c.f. ( p.l6| )). However, this is sufficient for our qualitative purposes of demonstrating 
the existence of constraints on the fermion mass implied by the boundary condition 



With these in mind, the boundary condition ( ^.131 ) reduces to: 



2 (3.14) 

2 Jo y^ + I [1 + mc^yy^ \ y^ + 1 J Sn^m^c V vr/ y 3 



To obtain the condition imposed on m by the boundary condition (|3.13| ) it suffices 



to observe that the first term on the right-hand-side is a function of mc^ alone, and 
that, after the (elementary) y integration, the resulting function of mc^ asymptotes 
rapidly to the value 3/4 (see figure |). 
This implies the following inequality: 

mG{0) < ^ fl - -) (3.15) 



37r V TT 

As already mentioned this bound is overestimated, given that in the actual situation 
the function G{y) is not exactly 1 immediately after the region y > m. 




Figure 4: Plot of the function f{z) = ^ Jq dy 
The function asymptotes rapidly to 3/4. 
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where z = mc^ . 



We next remark that m should actually be determined self-consistently from a 
solution of the pertinent gap equation. This will be done in section 5. However, 



at the moment, and for completeness, we shall assume that m is determined by its 
approximate form derived within the context of a large- treatment |l|] . Compat- 
ibility of the dynamical solution with the constraint ( ^.15| ) will then lead to further 
restrictions on the range of the allowed masses m. As we shall see in section 5, there 
is good agreement between the allowed fermion-mass ranges, obtained within the 
context of a large- treatment, and those obtained from a self-consistent solution 
of the mass gap equation within our approach. 

In the context of a large A^ treatment, and to leading order in 1/A^ resummation, 
the following solution for the dynamically-generated m is found f^, 

( . \ 



m 



C>(l)exp 



2n 



\ 



91 _ 1 

_7 J- 



(3.16) 



where g 
occurs [| 
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is the critical coupling, above which dynamical mass generation 
Compatibility of the solution ( p. 16 ) with the constraint ( p. 15 ) implies 
the existence of an upper bound on fermion masses, m < rrimax, where rrimax is 
defined through the intersection of the curves (|3.15|) and (|3.16|) in the (m, g) plane 
(see fig. This yields rrimax — 0.3. 
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Figure 5: Fermion mass versus the infrared-value of the coupling G(0). The solid 
curve represents the condition derived from the integral equation for the vertex, 
whereas the dashed line represents the solution obtained from the standard gap equa- 
tion in the large-N treatment. 



On the other hand, for large momenta, we know that G 1. Physically one 
expects a monotonic decrease of G{x) over the entire range of x G [0,oo). This 
would occur in our case if and only if G{0) > 1, which, in the context of the large- A^ 
result of 1^, |I|, implies a minimum bound for the fermion masses m > rrimin — 0.03. 
Actually, as we shall argue in the next section, G{0) should be comfortably larger 
than -y/3/2 for self-consistency of our approximations. 

Hence, we see that the monotonicity of the running coupling can be achieved in 
the context of a large- A^ treatment, if the mass m lies in the following regime: 

0.03<m<0.3 (3.17) 



or equivalently if the coupling at the infrared point G'(O) is restricted in the regime 
(see fig. §: 

1<G(0)<2.5, (3.18) 

At this point it is useful to turn to a numerical study of the equation (^.7|), supple- 
mented with the boundary conditions imposed by the solutions of the form ( p.ll|) , 
specifically: 



G(0) 




G'(0) 



8G3(0) 



(3.19) 
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where the constant c should be chosen in such a way that the bound ( |3.15 ) be 
satisfied. The numerical solution of ( p.7| ) is given in figure for a typical case, 
where m = 0.1, and G{0) = 2. As we observe, the figure clearly demonstrates the 
monotonic decrease of G in the small x interval, where equation ( |3.7| ) is valid. We 
also note that the solution asymptotes quickly to a constant positive value, smaller 
than 1 in contrast to ( ^.l(j| ), which would require G{x) ~ 1, for a; >> 1. This, 
however, presents no contradiction, because equation ( |3.7| ) was based on ignoring 
the inhomogeneous term Z{x). As one goes to a region of larger x, this assumption 
is no longer valid. In that case the inhomogeneous term Z{x) becomes important, 
and equation ( p.7|) receives additive (positive) contributions, resulting in G'(p) ~ 1, 
for large p, as required by ( |2.16|) . 
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Figure 6: Numerical solution of the equation ( \3. 1\ ) versus p/a, for a typical set of 

values m ~ 0.1 and G{0) = 2; the solution decreases monotonically and asymptotes 
quickly to a positive constant value. 



4 Equations for Photon and Fermion self-energies 

As a consistency check of our assumption about ignoring the effects of the inhomo- 
geneous term Z{p) as p — *• 0, next to G, we now turn our attention to the equations 
that determine the fermion and photon self-energies A{p) and B{p) respectively. As 
shown below, upon neglecting Z{p) in the infrared, these equations decouple from 



each other and both functions can be determined solely from knowledge of G{x). 
It should also be stressed that in the massless case the system does not admit a 
self- consistent solution. In contrast, the presence a fermion mass term changes the 
situation drastically by yielding self-consistent solutions for A{j)),B{j)) which are 
such that Z{p) 3-\/6/5, but at a rate slower than the one with which G ap- 
proaches G(0) as a; 0. Thus the approximation of ignoring the effects of Z{x) in 
the region a; << 1 is qualitatively correct. This is a highly non-trivial check of our 
approach, and justifies fully the approximations used above. 

To this end, we begin from the integral equation for which reads: 



6. 



which in terms of the semi-amputated vertex G becomes: 

d^k 



k^(k — py 



(4.1) 



(4.2) 



To arrive at the above equation we have carried out the 7-matrix algebra and we 
have used the approximation that G{k,p) = G{k). This approximation is equivalent 
to the one made for G in ( 2.12 ). It will be justified below, by the self-consistency of 
the solution. 

In terms of the dimensionless x, m parameters introduced previously, the equa- 
tion ([4.2|) becomes: 



X 



A-\x)x^-\ r dy 



y^G\y) 2x 



y2 _|_ 



T dy 



G\y) 

y2 _|_ yyip. 



(4.3) 



Differentiating twice with respect to x one obtains the following differential equation: 



d_ 

dx 



d 



X 



^(x'A-h 
dx ^ 



X 



8 G^jx) 







(4.4) 



Restricting our attention to the infrared region x — 0, one may ignore x^ in front of 
m^, and use the form ( |3.11D for G{x). After elementary integrations one then finds 
for A{x)~^: 



A- 
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-— X 
c' 
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— x"i + — x-^ln (1 



+ c?x 



1. , 3 



3 1 

-ex - -c^x^ + C(x^) 

Z: 



X 







(4.5) 



where c > is the integration constant of the solution ( p.ll| ), which obeys the 
restriction ( ^.15| ). We choose the constant of integration c' = 3/c, so that 



A{x^ 



,-1 



3 loo 

-ex C X + . . . 

5 2 



(4.6) 



as X — >■ 0. As we shall show later, this choice is compatible with the boundary 
behaviour Z{x) const, as x — >■ 0, advocated for the inhomogeneous term. 

Next we turn our attention to the equation for the photon function B{p), in the 
presence of a fermion mass. Following similar steps to the ones leading to ( [4.2|) , the 
pertinent integral equation reads: 



-'fJ.U 



P^S^u - Pf^Pu + B{p)e^ 



(27r) = 



G^ik)Tr 



{k — py 

(4.7) 

To trasnform the above equation into a scalar equation for B{p) it is necessary to 
take the trace on both sides. In doing so one assumes that the integral term on 
the right-hand-side is also transverse. This is expected from gauge invariance. In 
our truncated scheme it can be demonstrated that this is indeed the case, provided 
one makes our working hypothesis that G{k,p) = G{k). Specifically if one contracts 
both sides of ( |4.7| ) by p^, the right-hand-side vanishes (Ward identity for the photon) 
only upon making the above assumption for the semi-amputated vertex, and shifting 
appropriately the integration variable. 

After standard manipulations, similar to the ones above, one arrives at: 



d 

dx 



X 



dx \ dx ^ ^ y 
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G^x) = 0, 



X 







(4. 



Making the same approximations in the infrared x << 1 as in the case of A{x), and 
substituting the Emden's solution ( p.ll|) for G on the right-hand-side of ( [4.8|) , we 
obtain after some elementary integrations, upon setting their constants but one (ci) 
to zero: 

x^B~^{x) = -^x^ - 4^ + -^^^ + (l + c^a;) - ^In (l + c^x) + -cix'^ (4.9) 
c c c c c 8 



which for x << 1 yields: 



B-\x) ~ -|x"2 



1 X + — + - —x^ + 



(4.10) 



This form of B{x) gives rise to the following form for the photon propagator in the 
infrared /c — >• 0, 

A„.~^, k-^O (4.11) 



Notice, as a consistency check, that the photon continues to be massless in the 
chiral-symmetry broken phase. The corresponding static effective potential is given 
by the appropriate Fourier transform of the 00 (temporal) component of the photon 
propagator for ko = 0. In the case ( [4. 11] ) this yields formally an effective potential 
scaling like for large distances R, suggestive of confining behaviour The fact 

^However, this is a formal result, given that the corresponding momentum integrals are infrared 
divergent and hence need proper regularization, which falls beyond our scope here. For standard 
treatments see discussions in [Ol. 



that this behaviour is found in the chiral-symmetry broken case is consistent with 



general (four-dimensional) arguments [12, 14 that confinement is a sufficient but 



not necessary condition for chiral-symmetry breaking. 

An important feature of the expressions (^l6| ) and ( |4.8| ) is that, although they 
are derived only in case where a fermion mass m 7^ 0, however they do not explicitly 
depend on the magnitude of the mass. From ( |4.6| ) and ( [4.10| ) one obtains for Z{x) 
in the infrared region x ^ 0: 

Z{x) = B-'/\x)A-\x) ^ ^ - (1 + ^c^) + 0{x') . . . , a; ^ (4.12) 

where we have chosen the constant of integration ci such that there are no 0{x'^) 
terms. 

In order to check under which conditions the omission of the x-dependent parts 
of Z next to those of G{x) is justified, we must compare the corresponding linear 
terms in x of both quantities. Using that 

G(x) ~ (^m^c - m^c^x + . . . j , x (4.13) 

this requires: 



8/36-3(0) - 2G\0) - (37rV5)m2 » (4.14) 

This inequality furnishes a condition on G{0) under the assumption m << 1, which, 
as we shall see in section 5, turns out to be correct. The condition is 

G(0)»y|^1.22 (4.15) 

In the above formulas the symbol >> means actually at least an order of magnitude. 
As we shall see below, such a regime for G{0) arises self-consistently. Note that the 
condition ( [4.15| ) is in agreement with p.l6| ) given that Z{0) = 3^6/5 - 1.47 (c.f. 
Eq. (|4.12| )) is only slighty larger than ^3/2. 

In view of the existence of a non-zero constant value of Z{x) in the infrared 
region, our analysis on the restrictions (|3.15|) , (|3.18|) has to be repeated, given that 



these restrictions stem from the integral equation. From (^4.14|) , and the value Z{0) 
3\/6/5 it becomes clear that the boundary condition (|3.14|) is now modified to: 

12 mc^ /-i 1 1 / 2 \ 4/1 

+ ^ / dy—————— 1 + — — + — — - 1 - - 



Svrmc 2 Jo 7/^ + 1(1 + mc^y)'^ \ y'^ + 1 I Sir'^m'^c \ vr/ y 3 

(4.16) 

The first term continues to asymptote to 3/4, but now one obtains the following 
restrictions on the coupling: 

G(0) < 1^ + ! (1 - i) ^ (4.17) 

5 3 V vr/ irm 



From the above relation it becomes clear that, if G{0) < 12^/6/5 there is no restric- 
tion on m. However, for this regime of the couplings the condition ( 4.15|) , necessary 



for safely ignoring the effects of the inhomogeneous term Z{p) in the infrared next 
to G{p), is only marginally satisfied. On the other hand, if one allows G{0) to exceed 
this value, which is physically acceptable, then the following upper bound on m as 
a function of 6(0) is obtained: 



m < ^, G(0) > ~ 5.88 (4.18) 



G(0) 



5 



which replaces ( p.l5| ). 



One may repeat the comparison with the large N limit results 0, as in the 
previous section, to determine the new region of allowed values of the mass. The 
analysis is given in figure (|^). 
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Figure 7: Fermion mass versus the infrared-value of the coupling (5(0) for the case 
where the effects of the inhomogeneous term (almost constant) Z{x) have been taken 
into account in the infrared. The solid curve represents the condition derived from 
the integral equation for the vertex, whereas the dashed line represents the solution 
obtained from the standard gap equation in the large-N treatment. 



From fig. 1^ we then obtain the following allowed range of the coupling G'(O): 

5.88 < G(0) < 7 (4.19) 



to be compared with ( |3.18| ) associated with the monotonicity requirement for G{p). 



This analysis shows that the new allowed region for the fermion masses, in case 
( TO is: 

0.55 < m < 0.6 (4.20) 

where the lowest bound is the value of m corresponding to G = 5.88 from the large- 
A^-analysis formula ( ^.l(j| ). The region ( |4.2U| ) is to be compared with ( ^.17) ). As we 



shall see in section 5, the dynamically generated mass in our scenario, which does 
not resort to a large- iV treatment, lowers significantly the upper bound in ( |4.2U| ). 



It is important to emphasize that, as a consequence of ([4.5|),( ^rT0D , the non- 
amputated vertex Q{p), defined in ( p?7|) , approaches a non-zero finite constant value 
in the infrared p = 0. Notice that this constant behaviour is compatible with the 
Ward identity in the infrared limit p = 0, thereby implying that the information 
obtained about a non-trivial infrared fixed-point structure is gauge invariant, as it 
should be. 

This completes our analysis of including the effects of an almost constant Z{x) ~ 
Z(0) in the infrared. The above considerations justify retrospectively the omission 
of Z in deriving the differential equation (^3)- Indeed, a constant shift of G inside 
the derivative operators in that equation, which would account for a constant Z, 
does not affect its form. These results, therefore, constitute a highly non-trivial 
check of our approach. 



5 Dynamical Derivation of the Fermion Mass Gap 



In the previous section we have assumed the presence of a finite fermion mass, which 
we have treated effectively as an arbitrary parameter of the model. In this section 
we turn to the full problem, and study the dynamical generation of this mass, by 
deriving it self-consistently from the corresponding SD mass-gap equation. 

The equation for the gap is derived from the graphs of fig. 0, which yields: 



1 



Ji + M{k) ' (k-p) 



(5.1) 



where M{k) = Tj{k)/A{k) is the mass function, and we have pulled out factors of 
A{p) appropriately so as to be able to define an amputated vertex function G{k). 
The consistency of the approach will be demonstrated below by the existence of 
solutions. 

Taking the trace in the above equation, and performing the angular integration, 
one arrives easily at: 



M{p) = 3e^ 



Qa 
tt'^P 



dk 



(Pk 
(2^ 
k 



G\k) 



M{k) 



k'^ + M^{k) [k-pY 
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k^ + M^{k) 



G{kyM{k)\n 



k 



p 



(5.2) 



From the middle equation ( |5.2| ) it becomes clear that the quantity gii{k) = eG{k), 
defined in ( p.lO|) , plays indeed the role of a "running coupling". The situation is 
analogous, but not identical, to that of |TI|, where a "running coupling" has also 
been obtained from the gap equation in the context of a large analysis. However, 
in that case, one assumed an Ansatz for the vertex, satisfying a truncated form of the 
Ward identities. Instead, in our approach there was no necessity for a vertex Ansatz, 
since the non-perturbative vertex function was determined self-consistently from the 
SD equations. As a consequence, the running couphng ( |2.1CI| ) is constructed out of 



the amputated vertex function alone; the latter is a manifestly gauge-independent 
quantity, at least perturbatively. 

Passing into dimensionless variables, in units of a = M = M{k)/a, x = 

p/a, y = k/a, and working in the regime of low momenta x << 1 as usual, one 
obtains after some straightforward algebra, involving a truncated expansion of the 
logarithmic kernel: 



Mix) = 



x^Jo y^ + M^{y) Jx y^ + M^{y) 



(5.3) 



Differentiating twice with respect to x one arrives at: 



Given that dynamical mass generation is expected to be an infrared phenomenon, 
we now restrict ourselves to the region 0, 

x^ « « 1 (5.5) 

In this region we neglect next to in ( |5.4| ) and use the solution (|3.11| ) for 
G{x) ^ M^nc as X — > 0. The result is: 

(f ~ d ~ 

x—M{x) + 3—M{x) +9c^M{x) =0, x^O (5.6) 



Changing variables x — > ^ = x ^, the equation reduces to a Bessel equation |T7 

^'^^(0 - ^^^(0 + 9cV'M(0 = 0, e ^ oo (5.7) 
with (formally) general solution: 

M(0 = CieJ-2(-6cr'/') + C2iY^2{-Qci-"^) (5.8) 

where Ci^i = 1,2 are arbitrary constants, J_n(x) = (— l)"J„(x), n = 1,2,3,..., is 
a Bessel function of the first kind, and V„(x) is a generalized Bessel function of the 
second kind Jl^ . The latter is only defined for positive n and this imposes the choice 
6*2 = in ( [5.8[ ). Expressing the solution in terms of x, x — > 0, one has the following 
power series expression for the dynamical mass: 

oo 1 q 

M(x) = Cix-^ (3c)^"+^x^+" ^ c,c'+0{x), X (5.9) 

n=0 ^ \o -\- n) / 



^Note that this is the opposite Hmit than the one usuaUy considered in the framework of large- 
N analysis of dynamical mass generation, where one arrives at a gap equation by making the 
assumption a » p » m for the momenta p of the pertinent excitations. 



From this one obtains the following relation between G{0) and M(0) = mf/a: 



1 /3 

M(0) = rrif/a = (^^) Cl^^G^^'^O) (5.10) 



This result is to be compared the result ( p. 161) within the context of a large-iV 
analysis. In particular, at first sight it seems that the relation ( 5.10| ) does not 
have a critical coupling, above which dynamical mass generation occurs. However, 
because the result ( |5.10| ) has been derived in the context of the solution ( p.ll| ), 
one should bear in mind the restrictions characterizing that situation, in particular 
( [4.17| ). This implies an appropriate restriction for Ci 



This restriction implies a critical coupling, Gc = 12-\/6/5 ~ 5.88 but it is derived 
in a way independent of any large- analysis. The way to understand ( 5.11| ) is the 



following: one should first fix a range of G(0), with G{0) > 5.88, and then use a 
Ci that will be such that, within this range of the couplings, Eq. (|5.11| ) is satisfied 



for masses m « 1. As can be readily seen, the bound for Gi obtained from the 
requirement that m << 1 is far less restrictive than the one associated with ( |5.11 



provided G{0) is not too close to the critical G^ where the mass m vanishes. For 
instance, for G{0) = 0{8), the upper bound on Gi from (|5.11|) is of order C(10~^), 
while for G(0) = 6 the upper bound is Gi < 4. Notice that the bound is very 
sensitive to small changes in G{0). 

A typical situation is depicted in fig. ^ for two values of Gi = 10"^, 10~^. We 
observe that the case Gi = 10~^ yields an upper bound in the mass which is of order 
0.8 and hence should be discarded on the basis that it is not small enough. On the 
other hand, the value Gi = 10~^ yields an acceptable upper bound m ~ 0.1. In that 
case, from fig. ^, we observe that the allowed region of m is 

0.08<m<0.12 (5.12) 

to be compared with ( [4.20| ), derived using the result for the mass in the context of 
a large- iV analysis p. The corresponding regime of the couplings G{0) is: 

5.88 < G(0) < 11 (5.13) 

Before closing this section we would like to discuss possible applications of the 
above behaviour; specifically, the restriction ( [5.13| ) of the allowed values of G{0) 



may have interesting applications in case the above model turns out to describe the 
physics of high-temperature superconductors. We remind the reader that in such 
models dynamical mass generation coincides with superconductivity [Q-P]. In that 



The restriction ( ft. 17 ) is to be viewed as a boundary condition. We remind the reader that the 



requirement of finitcncss of M had already fixed the other constant C2 to zero. 
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Figure 8: Fermion mass versus the infrared-value of the coupling G{0) using l \5.1^ ) 
(dashed curves), for two values of Ci = 10"^ (lower dashed curve) and Ci = 10"^ 
(upper dashed curve). The continuous curve is ( U-lSj ), viewed as a boundary condi- 
tion. The value Ci = 10~^ should be excluded on grounds of yielding too high values 
of the mass m. 



case, the fermion fields ('electrons') of QED3 represent 'holons', i.e. electrically- 
charged excitations of fermionic statistics, which are constituents of the physical 
electron. The latter is believed to exhibit an effective spin-charge separation 
in the complicated ground-state of high-temperature superconductors. The 'pho- 
ton' of QED^ then represents an effective Heisenberg spin-spin antiferromagnetic 
interaction, responsible for binding the holons in Cooper-like pairs. In some mod- 
els 1^ the effective (gauge- invariant) coupling G(0) may be expressed in terms 
of the parameters of the microscopic condensed-matter lattice systems, whose long- 
wavelength limit is equivalent to the above QED^ model, as: 



J 



G(0)^~-(l-r/) (5.14) 



where 77 expresses the concentration of impurities in the system (doping), and J 
denotes the Heisenberg (antiferromagnetic) exchange energy. Hence, on account of 
( |5.13| ), Eq. ( |5.14| ) implies that 6 ^ (J/e^)(l— r^) ^ 11 for superconductivity to occur. 
In phenomenologically acceptable models [|| e^/J ~ 0.1, which implies an upper 
bound on ~ 0.4. However, the reader should bear in mind that the above-described 
limiting values are rather indicative at present, given that a complete quantitative 
understanding of the underlying dynamics of high-temperature superconductivity 
from an effective gauge-theory point of view is still lacking. 



6 Equation for the Vertex in the case of non-zero 
photon mass 



In this section we study the case where a (small) covariant photon mass 6 is added 
to the photon propagator ( |2.6|) : 

For the purposes of our analysis in this work, such a mass term will simply be added 
by hand, without further discussions about its origin. However, we point out for com- 
pleteness that a small photon mass may be the result of non-perturbative configura- 
tions (instantons) in compact U{1) three-dimensional electrodynamics |2^. More- 
over, a term acting like a photon mass term arises in real-time finite-temperature 
considerations in the latter case, of course, one loses manifest Lorentz covari- 
ance. 

It is straightforward to see that, in the presence of a photon mass, small compared 
to the scale a = the resulting integral equation for the amputated vertex reads: 



247r2 p 



(6.2) 



We are interested in the limit p << 6, which would allow us to study the effects of a 
photon mass on the infrared regime of the theory. Expanding the logarithm in (|6.2|) 
appropriately, we obtain: 



G(x) = — ^/ dyG\y) + — —^G^v) (6.3) 

where 6 = 6/a, x = p/a. 

Differentiating once with respect to x one obtains: 

^G{x) = -^G^{x)^^ (6.4) 

which can be integrated to yield 

G{x) = , ^ (6.5) 



c+ J^(^arctg(| 



X 



where c is an integration constant to be fixed by the boundary condition imposed 
by the integral equation (see below). 

The running coupling G{x) tends to G{0) = ^ a.s x —>■ 0. There is a non- 
trivial fixed point at x = 0, given that the renormalization-group /^-function ( |2.14| ) 
vanishes: 

d ^ . . 4 , , , , 



Notice that the function G{x) increases monotonically with increasing x, for small x. 
We next note that by choosing c appropriately one can satisfy the condition (|4.15|) 



G(0) >> y3/2, required for self-consistency of the approximation of ignoring the 
effects of the inhomogeneous term Z{p) in the infrared. In such a case one expects 
that at larger x the increase of the coupling stops at a certain (small) x = xq, 
and then the coupling starts decreasing to reach asymptotically the perturbative 
result G —>■ 1 (see fig. The existence of a local maximum is characteristic of 
the effects of the photon mass on the infrared behaviour, to be contrasted with the 
monotonically decreasing situation in the case of a non-zero fermion mass (see fig. 
^). We expect that in the general case, when both photon and fermion masses are 
present, the running coupling will exhibit a local maximum at low x when the photon 
mass is larger than the fermion mass, whilst this behaviour will be replaced by a 
monotonic decrease, in the case when the fermion mass is considerably larger than 
the photon one. This situation should be compared with the corresponding one in 



four-dimensional QCD in the presence of non- vanishing gluon and quark masses 114 
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Figure 9: Running coupling (versus momentum) in the case of a small photon mass. 
Note the hump at low momenta, which is absent in the case of non-zero fermionic 
masses. The dashed part of the curve is conjectural at present, and can only he 
derived numerically. Continuous curves are the result of analytic studies. 



Before closing we discuss briefly a physical situation which could be qualitatively 
similar to the case of a covariant infrared cut-off in the form of a non- zero photon 



mass, discussed in this section. This situation has been argued in ITT] to simulate 
finite-temperature effects, given that in such a case the photon propagator acquires 
a longitudinal plasmon mass term: 

A/,j.(po = 0, P 0, T) = /° + longitudinal parts (6.7) 



where, for simplicity, we restricted ourselves to the instantaneous approximation , 



Po = 0, and the plasmon mass term 



_ /2aln2T , , 

Moo(T) = y (6.8) 

We observe from the form ( |6.7D that the presence of the mass term Mqq behaves for 
low T somewhat analogously (but not identical to) a small photon cut-off mass S. 

Prompted by this observation we would like now to make some speculations 
regarding the low temperature effects of the S cut-off term in (|6.5| ). We hasten 
to emphasize that a discussion could at best only qualitatively correct, given that 
a quantitative understanding would require a finite-temperature extension of the 
above analysis. Due to the loss of Lorentz covariance in the respective propaga- 
tors the analysis is far more complicated than the zero-temperature one. However, 
when the temperature T is sufficiently low, its effects on the pole structure of the 
photon propagator may be simulated by a covariant photon-mass cut-off, at least 
qualitatively For this reason, we are motivated to examine the effects of 



5 << 1 on the low-temperature scaling of the resistivity |^ p; the latter is defined 
as the inverse of the conductivity o"/, which expresses the response of the system 
to a change in an externally-applied electromagnetic potential A^. Ohm's law, in a 



gauge with Aq = 0, then, gives |29, 11 



(P- + pl){l + Il{P.p„,T))^'--'' 

where Il{P,po,T) denotes the dimensionless photon polarization used in this work. 
Following the analysis of we may assume for our discussion below that [| : 



pl/a MlJa ~ 5^ ~ (6.10) 

TT 



From the discussion following (|2.10|) one immediately sees that the resistivity, defined 



for P = 0, = 5^, would be given in terms of the running coupling G[j)) by: 

p = 57(72(0, T) (6.11) 
Using ( |6.5|) in the limit = 5^, one has: 

The function p is plotted in figure (p!0|), versus a single-power scaling behaviour of 
the form: p ~ T^~^, with x small; the latter is characteristic of certain theoret- 
ical condensed-matter models [^, and has been used in order to bound possible 



^The analysis of is based on a real-time finite temperature approach. In such an approach 
the resistivity is defined in the infrared limit po ^ 0, P = in the presence of a plasmon mass 



term A/qo (6-8) in the photon propagator; as noted in ]2q |, the way in which the two mass scales 
{po,P) approach suffers from ambiguities, related to physical Landau-damping processes; it is 
not our purpose here to resolve such issues. 



deviations from the experimentally-observed linear-T behaviour of the resistivity of 
high-temperature superconductors at optimal doping [Q. 

We observe that, upon appropriate choice of the integration constant c, the 
temperature dependence of (|6.12| ) for a low-temperature regime, accessible to ex- 
periment, is hardly distinguishable from the single power scaling behaviour. 



0.8 



0.6 




Figure 10: Plot of the resistivity versus temperature T (in units of a), for a 

low-temperature region (continuous curve), and comparison with the experimentally 
observed single- scaling (basically linear) behaviour T^~^,x = 0.1 (dashed curve). 
The agreement is very good for a regime of temperatures accessible to experiment. 



7 Conclusions and Outlook 

In this work we have presented a study of chiral-symmetry breaking in QED^, 
based on a system of coupled non-linear SD equations. The novel ingredient is 
the introduction of a semi-amputated vertex, whose dynamics is governed by a 
suitably truncated SD equation of cubic order. This allows for a self-consistent 
determination of the infrared value of the effective charge, in the presence of infrared 
cutoffs provided by either the fermion mass or a covariant photon mass. The theory 
is characterized by a (gauge-invariant) non-trivial infrared fixed-point, suggestive of 
non- Fermi liquid behaviour |TT[]. 



The non-linear vertex equation furnishes highly non-trivial constraints among the 
infrared value of the effective coupling and the mass of the fermions. When these 
constraints are combined with the fermion mass-gap equation, they select specific 
regions of the coupling space for which dynamical mass generation occurs. It will 
be interesting to study how these constraints are affected if one goes beyond the 
one-loop dressed approximation considered here. 

From the physical point of view, we remark that dynamical fermion mass gener- 
ation in QED3 is associated with superconductivity fl], §]. In the relevant statistical 
models, whose low-energy effective theories are of QED3 type, the couplings depend 
on the doping concentration. Therefore the restrictions in the parameter space we 
have found above might impose restrictions in the allowed models. 



An interesting feature of our approach is the fact that the infrared structure of 
the photon propagator, derived as a consistent solution of a SD equation, imphed 
formally a confining effective potential. The fact that this behaviour occurs in the 
chiral-symmetry broken phase of the theory, appears to be in agreement with generic 
expectations that confinement is a sufficient condition for chiral symmetry breaking. 

An obvious next step in this programme is the attempt to solve the coupled 
system of SD equations numerically, without the approximation of ignoring the 
effects of the inhomogeneous term Z{p) in the infrared. We remind the reader 
that in the present work we have restricted the allowed values of the coupling such 
that this approximation was self-consistent. This, however, does not imply that 
solutions with Z{p) not meeting these requirements are impossible; to establish or 
exclude their existence one should solve the whole system of equations, which at 
present can only be done numerically. Such a numerical treatment will have the 
additional advantage of not resorting to the approximations made to the kernel of 
the integral equations in order to convert them into differential ones. 

Furthermore, an analysis at finite temperatures will reveal whether our conjec- 
ture in section 6, on the scaling behaviour of the resistivity, is correct. As we have 
discussed there, this issue acquires great importance in view of very accurate ex- 
periments in high-temperature superconductors indicating a basically linear scaling 
with temperatures of the resistivity in the normal phase (no fermion mass gap) of 
the system. 

Finally, an extension of these ideas to the non-Abelian case appears as a chal- 



lenge for the future, especially in view of recent claims [y, ^ that non-Abelian 



gauge symmetries might describe the dynamics of spin-charge separation in realistic 
antiferromagnetic condensed-matter systems p6| . 
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